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Abstract: We compute all iV-point primordial curvature correlation functions from inflation 
at tree-level up to N of order ten or more depending on the choice of parameters. This 
is achieved for resonant inflationary models in which the inflaton potential has a periodic 
modulation on top of a slow-roll flat term. These models find a natural UV completion in 
string theory implementation of axion monodromy. Key to the success of our computation 
is the observation that gravitational interactions among the perturbations can be neglected, 
which we argue is justified for any model of single-field inflation with parametrically large 
non-Gaussianity. We provide a comprehensive review and detailed derivations of known 
consistency relations for squeezed and collinear limits, and generalize them to any TV-point 
function. 
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1. Introduction 



Despite the increasing precision of recent cosmological observations, the data are still com- 
patible with Gaussian primordial perturbations as seeds of large scale structures. The most 
naive models of inflation, namely slow-roll inflation, generate perturbations that are Gaussian 
to very high accuracy. Nevertheless, Gaussianity is not a generic prediction of inflation, and 
indeed several models are known that lead to some observable non-Gaussian signal. A fu- 
ture detection of primordial non- Gaussianity would provide a powerful discriminant between 
different models. Given the wealth of new information hidden in the shape of the various 
iV-point correlation functions, non-Gaussianity also offers to precisely map the physics of 
inflation. 

For a Gaussian power spectrum, all iV-point connected correlation functions vanish. As 
a consequence, a simple way to measure non-Gaussianity is to look at higher point functions 
such as the bispectrum (three-points) or the trispectrum (four points). Given an inflationary 
model, it is often straightforward but computationally intensive to obtain these correlation 
functions. Most complications arise due to the coupling to gravity. For single field models of 
inflation, several bispectra have been computed (starting in earnest with [f, 2], see [3, 4] for 
reviews), together with a handful of trispectra [5, 6, 7, 8, 9, 10] and some partial results for 
N = 5,6 [11]. No higher iV-point function is known except in particular regions of momentum 
space 1 . 

The present paper stems from the following intuitive observation: when the most impor- 
tant interactions among the perturbations come from the sector driving inflation (typically a 
scalar field) as opposed to the gravity sector, then the latter can be neglected in computing 
correlation functions. More concretely one can formulate a "practical" decoupling limit in 
which the metric is an unperturbed (homogeneous) inflationary background 2 while the in- 
flaton is expanded as 4>(t) + 5(p(x,t). Any correlation function for 8<j) can be computed and 
translated into a correlation function for the curvature perturbations 1Z by performing a gauge 
transformation. This is very similar in spirit to [17], where it was noticed that this decoupling 
is analogous to the equivalence theorem in particle physics: W and Z scattering amplitudes 
at energies above their mass can be computed using only the Goldstone (longitudinal) modes 
[15]. A similar decoupling can also be obtained in the context of massive gravity [16]. Our 
observation is already implicit in the current literature and we will not attempt to formalize 

x See [12] and [13, 14] for some comments about the general structure of odd versus even iV-point functions 
in slow-roll inflationary models. 

2 Let us stress that this is different from usual choice of a flat gauge since in that case only the spatial 
components of the metric are unperturbed. 
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this decoupling limit. Instead, we exploit it to compute quantities that would otherwise be 
inapproachable. More concretely, we focus on resonant inflationary models and we find that 
upon neglecting gravitational interactions the calculation simplifies to the point that we are 
able to compute all iV-point curvature correlation functions at tree level for N as large as ten 
or twenty depending on the parameters. 

Resonant models were first proposed in [18], where it was noticed that in the presence 
of a periodic modulation in the inflationary potential, non-Gaussianity can be resonantly 
enhanced. These models found a natural UV completion 3 in string theory constructions 
based on axion monodromy [20, 21, 22, 23]. The basic idea is to consider an axion that enjoys 
a shift symmetry to all orders in perturbation theory and use it as the inflaton [24]. This 
symmetry can be broken by non-perturbative effects but also explicitly by some small effect. 
The latter possibility provides a mechanism to naturally protect the flatness of the inflaton 
potential in large field models of inflation. Since non-perturbative effects are generically 
present, one expects a potential of the form 



where is a canonically normalized axion with axion decay constant /, A is some non- 
perturbatively generated scale and V sr is an arbitrary slow-roll potential proportional to the 
explicit breaking of the shift symmetry. A potential of the form (1.1) has been obtained in 
some class of type IIB string theory compactifications [21, 22, 23], where V sr ~ ficj). In the 
following we will not assume this specific linear dependence, but consider the more general case 
where (1.1) is monotonic, which translates into the inequality b < 1 if we define b = A 4 /(/V^,' r ). 
The classical evolution induced by (1.1) has two typical time scales: the Hubble parameter 
H and u= \4>\/f whose ratio we denote by a = uj/H. 

In the regime a> 1, non-Gaussianity is greatly enhanced by a resonance between the 
perturbations and the background [18, 25, 26]. Since non-Gaussianity is parametrically large 
in a, the decoupling of gravity can be used and this drastically simplifies the calculation of 
correlation functions. In fact, we are able to compute the single- vertex contributions to any 
(connected) curvature correlation function (iZ V i.e. for any N, at tree level and away from 



3 For another occurrence of oscillations in string theory models see [19]. 
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collinear singularities. Working in an expansion in small b and large a, we find 

(K kl ■ . .^ kjv) -Sle vertex = ^3^3 A N B N (h) , (1.2) 

An = ( _)^:!^ a 2iv-9/2 A ^-2 (2 ^ J V-l 



f 



where K = ^ = 4>* — \/2e* In K/k*, k* is some pivot scale and <^>* the corresponding 

value of the scalar field. The power spectrum is well approximated by the slow roll result 
A|j = with e the slow-roll parameter. Notice that the shape only depends on the 
magnitude of the vectors but not on their respective orientation. This gives a total of N 
variables as opposed to the usual 3N — 6. Note also that the above result reproduces the 
bispectrum calculated in [26] using the full perturbed metric. 

We provide convincing evidence that the single-vertex diagram gives the largest contri- 
bution to any iV-point function away from collinear limits up to some iV max which depends 
on the parameters of the model (as summarized in figure 3), i.e. 

<rc kl • ■• Wk W ) M v * ■ ■ .^ k j* ertex (L3) 

JV <-/V ma x 

As figure of merit, for a = 100 and any b < .1, N max = 9, while for a = 1000 and any b < .1, 
Af max = 17. For iV > N max the number of multi- vertex diagrams is so large that collec- 
tively they become comparable with single- vertex diagrams, although they are individually 
much smaller. Computing all these diagrams is straightforward in principle, but becomes 
computationally intractable for large TV. 

One potential concern is the behavior of (1.3) in specific momentum configurations. 
Generically, an iV-point function has singularities when one of the momenta (squeezed limit) 
or a combination thereof (collinear limit) is taken to zero. The squeezed limit is correctly 
captured by the subleading term in (1.2) but collinear limits are not reproduced. Fortunately, 
relations are available to obtain the behavior of correlation functions near these singularities as 
first pointed out in [2] and subsequently studied and generalized in [27, 6, 53, 8, 29, 30, 31, 32]. 
We review the known consistency relations and generalize them to any N. We provide a de- 
tailed discussion of their derivation, both in momentum and coordinate space in appendix B. 
Using these relations we can supplement (1.3) with the leading divergent terms near collinear 
singularities. As an example we explicitly write down the resonant trispectrum in section 6. 
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The outline of this paper is as follows. In §2 we explain our general strategy of decoupling 
gravity. We then apply it to the resonant inflationary model which we introduce in §3. In 
§4, we compute all single-vertex diagrams. We show that the first two leading terms in a 
are compatible with known results and that the consistency relation in the squeezed limit is 
obeyed for any N. We then estimate the size of multi-vertex diagrams by using consistency 
relations in the collinear limit in §5. We combine all our results together and discuss their 
limits of validity in §6. We conclude in §7. The appendices provide the Feynman rules for 
the theory under consideration and review and prove various useful consistency relations. 

2. The General Strategy of Decoupling Gravity 

In this section we outline the general strategy we will adopt in this paper for the computation 
of correlators of primordial perturbations. We are surely not the first to realize and use this 
strategy but we would like to make the statements that are often implicit in the literature as 
clear as possible. 

A typical theory of inflation may look like 4 



with some Lagrangian describing the sector responsible for the energy momentum tensor 
that drives inflation, e.g. a scalar field with a slow-roll flat potential. The standard procedure 
is to first find a quasi de Sitter solution of the classical equations of motions derived from S. 
One then expands the action in powers of the perturbations of the inflaton and the metric. In 
the presence of gravity this is a laborious computation due to the large number of interaction 
terms. 

As we mentioned in the introduction, we observe that 

• Single field models of inflation with parametrically large non-Gaussianity admit a limit 
where the leading contribution to curvature correlation functions can be calculated using 
only the action for 5(j) while leaving the metric unperturbed. 

By parametrically large, what we have in mind is that in the given model parameters exist 
that we can tune to make higher Appoint functions arbitrarily large compared to the power 
spectrum which is held fixed. The basic idea is that gravitational interaction are naturally 
small in the context of inflation and therefore if there exists an independent parameter which 

4 We set the reduced Planck mass Mp1 d to one. 
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allow large interactions it should be possible to neglect gravity. The end result is expected to 
be correct only to leading order in that parameter. 
The computational prescription is the following: 

• Solve the Friedmann and scalar field equations to find the classical background, say 
{a(t)^(t)}. 

• Derive the interaction Hamiltonian Hj for the scalar field perturbations 8(j) from the 
action S = J ^—gC^, where g is computed from the background, e.g. \f--g = a(t) 3 . 

• Compute the equal-time correlators (5(j)(t) N } in the in-in formalism using Hj computed 
as above. 

• Convert the correlators of scalar field perturbations computed at horizon exit into cor- 
relators of comoving curvature perturbations 1Z, which are frozen outside the horizon. 
This conversion can be done at linear order in perturbations, i.e. using 5(j) = \j2e^R, 
(where e* = —H/H 2 is the Hubble slow-roll parameter). Note that to ensure that the 
gauge transformation is linear one may have to perform it soon after horizon crossing. 

The resulting correlators (1Z ) should be accurate at leading order in the parameter which 
makes non-Gaussianity large. 

Before we plead further, we would like to make the following disclaimer. The idea that 
correlators of perturbations might be dominated by interactions inherent to the inflationary 
sector and that gravity might give subleading corrections is certainly not new. Well known 
examples are given by models with small sound speed and/or derivative interactions such 
as DBI inflation [33], p-adic models [34], ghost inflation [35], quasi-single field models [36] 
and galileon[37]. In most cases, the authors computed the non-Gaussian signal following the 
prescription outlined above (see for example [38, 39, 40] and [6, 7, 9, 17]). 

The purpose of the present paper is not to give a general proof of the validity of the above 
strategy, but to use it to compute quantities, such as iV-point primordial correlation func- 
tions for N 3> 4, that are yet unknown and unthinkable to compute including gravitational 
interactions. The success of our computation provides additional motivations to both put the 
decoupling strategy on firmer ground, which we hope to come back to in the future, and to 
further employ the decoupling strategy to widen our knowledge of correlation functions from 
inflation. 
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3. Background and Power Spectrum 



In this section we start up with a review of the model and of the classical inflationary evo- 
lution that it generates. This constitutes the background for the computation of correlation 
functions. We then derive the spectrum and show that it agrees with previous results [22, 26]. 
For more details about the phenomenology of the spectrum and the bispectrum we refer the 
interested reader to [21, 22]. 

Consider a modulated potential of the following form 

Vtf) = y sr (<A) + A 4 cos(0//), (3.1) 

where V sr {4>) is a smooth potential that, by itself, would admit at least 60 efoldings of slow-roll 
inflation and A and / are two parameters with dimensions of mass. The sinusoidal term arises 
naturally in the presence of non-perturbative corrections to the inflaton potential, in which 
case A is some non-perturbatively generated scale. A familiar example in field theory is an 
axion that couples to some gauge sector with strength given by the axion decay constant / 
5 . Instanton configurations in the gauge sector can induce periodic corrections to the axion 
potential, which have periodicity 2tt/ for a canonically normalized field and A oc e -1 / ff2 with 
g the gauge coupling. An axion then is a good inflaton candidate because it enjoys a shift 
symmetry to all orders in perturbation theory. Even when the shift symmetry is explicitly 
broken in a controlled way, e.g. by the V^ r term in (3.1), its remnants still protect the flatness 
of the potential against quantum corrections. This provides a technically natural realization 
of large field inflation. It has been argued [21, 22] (see also [43]) that a certain class of string 
theory constructions leads to potentials of the form (3.1). For the purpose of this paper we 
take the potential (3.1) as given and study its phenomenological consequences. 

With a suitable V ST , which we assume, the potential (3.1) can drive a prolonged phase of 
exponential expansion of the universe. The classical homogeneous background resulting from 
solving the Einstein and scalar field equations at leading non-trivial order in A was computed 
in [22]. It turns out to be useful to introduce the monotonicity parameter 

For b < 1 the potential (3.1) is monotonic (provided V sr is, which we are assuming). Com- 
parison with observational data 6 , e.g. with the power spectrum of the CMB, tells us that 



5 The consequences of this coupling for inflationary perturbations have been studied in [41] and more recently 
in [42]. 

6 Strictly speaking this comparison have been carried out only for two explicit forms of V BI , linear [22] and 
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b <C 1. This justifies working in an expansion in powers of b. Another restriction on our 
parameters is / <C Mpi. This choice is motivated by the fact that, despite many efforts, su- 
perPlanckian axion decay constants have not been realized in any string theory construction 
so far and there are general arguments [46] that they can not arise in controlled situations, 
i.e. within the supergravity approximation with a perturbative string coupling. When b <ti 1 
and / <C Mpi, the sinusoidal oscillations in the potential (3.1) induce a time dependence of 
the background in addition to the one due to slow-roll inflation. We call oj the time frequency 
of the oscillations of the background, and we introduce the dimensionless parameter 

where e* = —H/H 2 is the first Hubble slow-roll parameter (evaluated at some pivot scale). 
Throughout this paper we will assume to be in the regime a>l, i.e. that the time scale of the 
oscillations of the background is much shorter than the time scale of the slow-roll inflationary 
evolution. This is the regime in which cosmological perturbations are resonantly enhanced, 
as first noticed in [18], by the interaction with the background. The study of this effect is the 
main focus of the present work. 

The solution for the background at leading order in b and for a ^> 1 is [22, 26] 

«i)-A,M-3g *.(*#>), (3.4) 

where (po(t) is the time evolution of the scalar field <j) without the sinusoidal modulation, 
which of course depends on the specific choice of V SI (4>). For our purposes we do not need to 
specify V sr (4>), so our results are valid for any slow-roll flat potential V SI - In the following it 
will be useful to rewrite (3.4) as 

<p(X) = 4>k + y2e* In A y= sin 



2e* V / 

4>K = 4>* - y/2e^kiK/k* , (3.5) 

where time is now parameterized in terms of X = —Kt, with adr = dt being the conformal 
time, K is arbitrary (it cancels out) and is some pivot scale relevant for the experiments 
under consideration. 

quadratic [44]. On the other hand we expect qualitatively similar results for other slow-roll potentials V ST - For 
a partial list of constraints on oscillations in the CMB see [45]. 
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3.1 The Two-point Function 

In this section we compute the two-point function of primordial curvature perturbations 
following the general strategy outlined in section 2, i.e. using only the scalar field sector 
of the theory. We show that at leading order in b <C 1 and a 3> 1 this agrees with the 
computation performed in [22, 26], where gravity was taken into account. 

Before proceeding, let us review the definitions of the degrees of freedom we will work 
with and our gauge choices. Generic scalar metric fluctuations, to first order in perturbation 
theory, can be written as 

ds 2 = -(1 + 2§)dt 2 + 2a{t) 2 B )i dx i dt + a(t) 2 [(1 - 2¥)<5y + 2E#] dx i dx j , (3.6) 

where a(t) is the scale factor and {<!>, B, ^, E} are scalar perturbations. At first order in 
the perturbations, after imposing the constraint equations and in the absence of anisotropic 
stress, the five perturbations (four in (3.6) plus 5 (ft) reduce to just two, that combine together 
into 

Tl = -^f - ^S(ft , (3.7) 
(ft 

which is invariant under time diffeomorphisms. In comoving gauge the only degree of freedom 
is in the comoving curvature perturbations 

6(ft = , 9ij = a 2 e 27 % . (3.8) 

With an appropriate time diffeomorphism one can go from comoving (or uniform density) 
gauge into spatially flat gauge where = and vice versa. Since 1Z is gauge invariant the 
two gauges are related by y/2e*ty = — 8(ft where y/2e^ = \(ft\/H. 7 

Let us now consider the linearized classical equation of motion for the scalar field pertur- 
bations in momentum space S(ftk (due to rotational invariance of the background 5 (ft depends 
only on the absolute value k of k) in spatially flat gauge: 

S(ft k + K,5(ft k + 3H6(ft k - 3bH 2 a5(ft k cos ( ^f] = , (3.9) 



a \ J 

where the slow-roll part of the potential has been neglected. Going to comoving gauge, after 
a long but straightforward calculation one can rewrite the equation above as 



d 2 x n k - 2{l + 5 + e * ] d x n k + n k 



. , d 2 V2^ 3b 

1 H — ^QCOS . 

x 2 V / 



0, (3.10) 



7 In this paper we assume that (j> < 0. If the opposite is true, the sign of the gauge transformation will 
differ. 
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where x = —kr and 5 is a Hubble slow-roll parameter 8 denned by 5 = H /(2HH). Because 
of the cosine part of the potential we have the following hierarchy of slow-roll parameters 
e* <C 5 <C 5/H. Now notice that 



1 



H 



+ 5(l + <5 + 3e*) + e*(l + 2e*) 



x 2 H 



(3.11) 



In the background induced by the potential (3.1), one finds [22] 



■ ~ 36a cos 
H \ J 

Hence (3.9) agrees with eq (2.17) in [26] 



Tixx ~ 2 ^ 1 + ^ K X + K = 0. 

x 



(3.12) 



(3.13) 



at leading non-trivial order, i.e. for e* <C 5 <C 5/H. 

To summarize, we have substituted the relation 5(f) = y/2e*lZ into the equation of motion 
for (f> and showed that the resulting equation is the same as the one used in [26, 22] at leading 
order. Since the power spectrum follows directly from the solution of the Mukhanov- Sasaki 
equation (3.13), we conclude that using just the ^-sector and neglecting the interactions 
coming from gravity is enough to correctly reproduce the power spectrum with its oscillatory 
features, to leading order. 

The quantization of the perturbations proceeds as usual and we briefly sketch it in order 
to introduce our notation. The momentum space perturbations are related to the real space 
perturbations by 



<ty(x,t) 



d 3 k 



5<f>(k,t)e 



ikx 



(3.14) 



where the three- vectors {x, k} are the comoving coordinates and momentum, respectively. 
5(j)(k, t) is determined by the solution 5(f>k{t) of (3.9) with Bunch-Davies initial conditions. 
Then the momentum space quantum field is given by 

<ty(k,t) = <fyfc(t)a k + tyfc(t)aL k , 
where the creation and annihilation operators satisfy the commutation relation 



Ok,aL k . 



(2vr) 3 <5 3 (k + k / ) 



A summary of various definitions of slow-roll parameters and their relations can be found in the appendix 
of [26]. 
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Parameterizing time with x = 



—kr one finds 



6(j) k (x) 



H 



-iH 



2 3/2 

l-ta:)e fa + 0(e.,<5, 6) 



(3.15) 



Note that the mode functions <5</>fc are given by the standard de Sitter result plus two types 
of subleading corrections. First, slow-roll suppressed terms that would change the Hankel 
function into H„ with v = 3/2 + 2e* + 5. These terms would correct our results for non- 
Gaussianity by a small logarithmic running so we neglect them from now on. If needed it is 
straightforward to keep these terms. Second, b suppressed terms coming from the oscillation 
of the background. These terms are important to compute the leading correction to the 
slow-roll power spectrum, which turns out to be 



\R k \ 2 



(27r) 3 ^ 3 (k + k')|^| 2 



2ir 2 



n 
k 3 



1/2 



0k 



1 - :',!> ( ^ J cos ^ f 



(3.16) 



where A J? 



H 2 



v . T 2 e - On the other hand, as discussed in [22, 26], for the bispectrum and higher- 
point correlation functions these terms can be neglected. In fact, the leading term at linear 
order in b in the correlators arises when taking the 6° terms from the mode functions and the 
b term from Hj (which of course does not appear in the computation of the spectrum) for 
the mode functions. 



4. Single-vertex Diagrams 

The iV-point correlation functions of (50k can be computed using the in-in formalism. We will 
perform the calculation using the canonical formalism, but we will organize it according to 
the structure of (in-in) Feynman diagrams following the notation and Feynman rules in [47]. 
In this section we compute the contribution to iV-point scalar correlation functions coming 
from diagrams that have a single interaction vertex or equivalently no internal propagators. 
As an example, the corresponding Feynman diagrams for N = 3, 4 are depicted in figure 1 
(the relevant Feynman rules are collected in Appendix A). Multi- vertex diagrams will be 
discussed in the next section. We follow the strategy outlined in section 2, i.e. we only use 
the scalar field Lagrangian. 
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For our canonically normalized field 5(f) the interactions comes purely from the potential, 
hence the Hamiltonian is 



/oo 
d 3 xa' 

N=3 



6<p(x 



(4.1) 



Since we are interested in tree-level single-vertex diagrams we start from the leading order 
result in perturbation theory 



II^(r)\ = if_^ ^ <[50 kl • • ■S^ n (t),H i (t% 



(4.2) 



Here we have assumed a quasi de Sitter background where, to leading order in the slow-roll 
parameters 9 , the scale factor is given by a(r) = — Note that in order to simplify the 



+ 



x 




Figure 1: Single- vertex Feynman diagrams for three- and four-point functions. See Appendix A for 
the Feynman rules. At N = 5, there are three diagrams and so on. 

notation, in (4.2) we have denoted the connected part of the N-point function by () instead of 

9 We are allowed to use the slow-roll approximation for the background even in the presence of oscillations 
because we are working at leading order in b. The largest contribution to the correlator comes from taking 
the linear in b term in Hi while the b° terms from all other factors. 
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() . In this paper we only refer to connected correlation functions so we always drop the label 
since no confusion can arise. To compute the commutator in (4.2), we expand the real-space 
fields 5(p(x) appearing in Hj in Fourier modes. Since it is 7Z and not 6(p that is constant 
outside the horizon, it is customary to evaluate (5cj)(T) N } at the time of horizon crossing and 
perform a gauge transformation at that time to the variable 1Z. This is not crucial since the 
evolution of 5<p is a small slow-roll suppressed effect, which should precisely cancels the time 
evolution of the factor ^/2e^ involved in the change of gauge. As compared to the effect of 
oscillations, these effects together with other slow-roll suppressed corrections are subleading 
and we neglect them throughout our computation. Hence, it is consistent to focus on the 
limit t — > where the relevant commutator is 



N N 

n^(°)>n^( T ') 

1=1 j=i 



tt2N 



N 

AKt 



J(l - iferO^Cki - ki) - c.c 



i=i 



(4.3) 

where K = Y2i ki an d we have subtracted the divergences by normal ordering. Putting this 
back in (4.2) and performing the d?ki and d?x integrals gives 

< 1P<M0) ) = =-r p In, (4.4) 



>° « 

in 



^-{-i)VW [(i _ i klT ') • • • (1 - ik N r')e iKT ' - c.cl . (4.5) 



Now normally V^ N ' is a constant and can be taken out of the integral in which cases one can 
show that [48] 

"YklNe, (4.6) 



IN 



2 



y(JV) 3 

where N e = ln(— Kr) is the number of efoldings. For slow-roll models, the contribution from 
this part of the Hamiltonian is actually slow-roll suppressed compared to the coupling to 
gravity so this is usually dropped (the iV e enhancement disappears when converted to 1Z in 
single field models of inflation). For resonant models, this term gives a large (often leading) 
contribution to iV-point correlation functions as we will show. 

In order to compute the integral (4.4) we have to specify V^ N \ Let us consider the 
potential in (3.1). The slow-roll 10 term V SI gives rise to interactions much weaker than those 
coming from the the oscillatory part and can hence be neglected. Therefore, to leading order 



10 By "slow-roll" we mean not only the condition e,j< 1, but also that all higher slow-roll parameters (e.g. 
see [49]) are much smaller than one. In the presence of sinusoidal oscillations the first condition is satisfied 
while the second is not. 
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in the slow-roll approximation, we consider 

V{4>) = 3a6/ 2 F 2 cos(0//), (4.7) 

where the monotonicity parameter b was defined in (3.2) and a in (3.3). Expanding this 
potential in perturbations we find the following iV-th derivatives 

^^fH^cosW/fl.^even 
f \(-)^sin(0//) , TV odd, 

In the above expression the time dependence is in H and <j). On the inflationary background 
we can use 4>(X)/f = 4>k/ f + alogX, where we parametrize time with X = —Kr and K is 
arbitrary at this point. After expanding out the products in (4.5), we find integrals of the 
schematic form 

J dXX n sin (X + 9 X ) sin (a log X + 9 2 ) , (4.9) 

for some integer n and some real phases 9\ and 9 2 • The phase of the integral is stationary at 
X = a and this is where the integral gets its main contribution. These kind of integrals can 
be performed using the stationary phase approximation 11 

/oo 
G(x)e ikf( - x) dx = G(x)e lkf W 
-00 



(4.10) 



where x is the point where the phase f(x) becomes stationary, i.e. f{x)' = 0. For example 
using this formula one can compute 



// 2"7T 
dXf(X) cos(a + a InX) sin(X) f(X)^-a 1/2 sm(a - 6) , 

dXf(X) sin(a + a In X) cos(X) ~ f(X)^-a 1/2 sm(a - 9) , 

dXf(X) cos(a + a In X) cos(X) ~ f(X)^-a 1/2 cos(a - 9) , 

dX/(X)sin(a + alnX)sin(X) ~ f(X)^-a 1/2 cos(a - 6») , (4.11) 

where the phase is stationary at X = X = a and 9 is a fixed (irrelevant) phase. We will 
have a sizable effect if a S> 1, which physically corresponds to the regime when the time 
scale of the oscillation of the background is much shorter than the time scale of the slow-roll 
evolution. In this regime, the leading term is the one with the highest power of X. 

11 In fact, a direct integration is possible as well and gives some combinations of F functions. On the other 
hand, at the order we are working these V functions should be expanded using Stirling's formula such that the 
end result is the same as the one from the stationary phase approximation. 
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4.1 The Leading Contribution 

The whole integral (4.4) can be computed (leading to (4.20)), but here we proceed step by 
step and first compute the leading piece. The leading term in powers of a comes from the 
highest power of X in the integrand of (4.4) 

I N = ~i dX V^ N \X)X N - i g{k i ) [[(i) N e~ iX - c.c] + . . . , (4.12) 
Jo 

where the dots stand for terms suppressed by at least one power oil/a and we have introduced 
the notation 

Q(h) = ^. (4.13) 
Plugging the explicit expression for (4.8) into In as given in (4.12), we obtain 

3abH 2 , s f°° at-4 J cos ( + oilnX ) , N even 



lN = -2-^g{h) dXX^-U } J { . (4.14) 

f N 2 Jo I cosXsinf ^ +alnX) , iV odd 

Using (4.11), we find that both integrals give the same answer up to a sign. Modulo a 
^-independent irrelevant phase, the result is 

In = {-) N V^j^G{h)a N ^ sin . (4.15) 

We can now write the leading contribution to the curvature correlation functions by perform- 
ing the gauge transformation ^2e*lZ = 5(f). We find 

where 

A N ^ (-) N3 ^a* N ^(2^Al) N ^ , (4.17) 

with = g^-- It is straightforward to check that for N = 3 this result agrees at leading 
order with the one obtained in [26] , where also the interactions induced by gravity were taken 
into account. 

4.2 Main Result and Subleading Contributions 

It is also interesting to compute the first subleading terms, i.e. those with the second largest 
power of t' (or equivalently of X) in (4.4) since they are formally divergent in some squeezed 
limits. Like for the leading term, these terms agree with [26] for N = 3. The computation is 
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completely analogous to that of the leading contribution, so we give directly the final result 
of leading plus first subleading term: 
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(4.18) 



where F{ki) = {jXj^i ^jj K^~ N and we have omitted a /c-independent phase in the 
argument of the trigonometric functions since this is arbitrary in the model. Putting things 
together we find 
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(4.19) 



This is our main result. It is natural to compare this result with the one in [26] at subleading 
order as well. The two look very similar except that the sum in (4.19) runs over all i and 
j, while in [26], e.g. in equation (1.3), the corresponding sum runs only over i ^ j. The 
difference is given by — 3 cos {4>k/ f) /ot. This term is nothing but a different choice of the 
fc-independent phase omitted in the trigonometric functions in (4.19). In order to see this, 
consider shifting <Pk/ f i n (4.19) by the ^-independent phase 3 /a. By using trigonometric 
relations and expanding in ft ^ 1 one finds exactly the result of [26]. Since the phase is 
anyways arbitrary, we conclude that (4.19) exactly reproduces the results of [26] both at 
leading and at first subleading order in a. 

It is straightforward to compute all other terms in (4.4). We will use these terms in 
section 6 to estimate the combinatoric coefficients of different diagrams. The result is 

l 



3abH 2 a N ~ 4 A 

In = (-) V2Tra jAf _ 9 „ N _z l^ 
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where I = is the leading term. When we put this back into (4.4) we find 
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(4.20) 
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As we just discussed, the I = 0, 1 terms exactly reproduce the result of [26]. On the other 
hand, at the next order in the large a expansion the calculations differ. This is presumably 
due to the gravitational interactions, which we have neglected and which contribute to the 
terms I > 2. Despite this, (4.20) is still very useful to estimate the combinatorics factors 
arising for large iV as we will explain in section 6. 
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4.3 Squeezed Limit 

Interesting limits and consistency relations arise when some set of momenta goes to zero. 
We will be interested in two classes of consistency relations: in this subsection we consider 
relations arising in the squeezed limit, i.e. when one of the external momenta goes to zero. 
In the next section we will make repeated use of a relationship arising in the collinear limit, 
i.e. when a (proper) subset of momenta sum up to zero. We present an extensive review of 
these consistency relations and how to derive them in appendix B. In the case of the squeezed 
limit, we can perform a non-trivial check of (4.19). 

In the limit k\ <C &2 ~ &3 • • • ~ k]y, the Appoint correlation function should obey 

a 



Urn (K kl 

fci->0 



K kN ) = (2ir) 3 5 3 (K)\K kl \W2e*--(K ]i2 ---K k! 



(4.21) 



where (} indicates that the correlator does not contain the factor (2-7r) 3 5 3 (K). Although this 
is a relation between the full N- and (N — l)-point correlation functions, it should be satisfied 
order by order in an expansion in b and a. We will use the leading term in b and a from 
(4.19) in both the right and left hand side of (4.21). Let us start from 



(n N )' = a n 



K 3 ' 
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sm(a) cos( 
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(4.22) 



where we remind the reader that a = 0k/ f- Then, to leading order in k\,b and a, (4.21) 
leads to 



A N 1 

A N . X Kkf 

or simply 



1 I \ L-ii hi 
■— cos(a) — ^ — 
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— - | = 2vr 2 ^acos(a) , (4.23) 
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(4.24) 



which is an identity given A^ in (4.17). Hence we have shown that the relation (4.21) is 
satisfied to leading order in k\,b and a for any N. 



5. Multi- vertex Diagrams 

In the previous section we have computed all tree-level single-vertex diagrams at leading 
order in small b using the strategy outlined in section 2. The final result is given in (4.20). 
Here we discuss the much harder problem of multi-vertex diagrams. We will not directly 
compute them, but we will estimate their size using some relations arising in the collinear 
limit. Details of the derivation of these relations can be found in appendix B. The results 
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collected in this section provide us with an expectation of when multi-vertex diagrams are 
important as compared to single-vertex diagrams. On general ground, multi-vertex diagrams 
are suppressed by extra factors of b and a but their contribution can still be important, if 
there is a large number of terms and/or around singular points. We will find that for the 
first few correlation functions (small N) multi-vertex diagrams are subleading, while they 
grow fast for larger N. Combining the discussion of the current and previous sections into an 
accurate estimate of TV-point correlation functions will be the subject of section 6. 




Figure 2: All tree-level two-vertex diagrams for the trispectrum up to permutations. These diagrams 
are subleading to the one computed in figure 1 for small b and large a. In the collinear limit where 
q = + k 2 — > these diagrams diverge as q 3 and will eventually dominate over the single- vertex one 
in (4.19). 



5.1 Collinear Limit with Scalar Exchange 

The collinear limit is the limit in which a (proper) subset of the momenta {ki, k2, . . . , k/} for 
1 < I < N of some Appoint correlation function sum up to zero. In this limit, the Appoint 
correlation function is dominated by diagrams describing the exchange of a long-wavelength 
(momentum q — > 0) classical scalar and tensor mode between an I- and an (N — Z)-point 
correlation functions as depicted in figure 5. Here we start by analyzing the scalar exchange. 
In Appendix B we find 

where q = Yl\=i denotes "scalar exchange", and we remind the reader that the prime 

in ()' means that the factor (27r) 3 5 3 (K) is omitted. The relation (5.1) should be verified order 
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by order in b. Since we have computed the single-vertex diagrams for any N we can use 
them in (5.1) to estimate the two- vertex iV-point function in the collinear limit in which the 
propagator between the two vertices goes on shell. 

As an example let us consider N = 4, i.e. the trispectrum, which contains three tree-level 
multi- vertex diagrams depicted in figure (2). If we estimate them in the collinear limit by 
mean of (5.1) we find 

• " n ^ fE = 2e * |7 ^ |2 |; <^>' h {nk)l ' (5 ' 2) 

where fcy = |ki + kj|. Using the power spectrum to leading order in b given by (3.16), the 
trispectrum in the collinear limit becomes 12 

< k - ■ ■ ■ n "- fE = ^ si " 2 a <2 2^r • <5 ' 3) 

Two comments are in order. First, we notice that although the trispectrum depends in general 
on six real parameters (i.e. 3 x 4 = 12 momentum components minus three for translational 
and three for rotational invariance), in the collinear limit, to leading order in k\2, it depends 
only on three, namely ki2,ki,k%. This is because, since the exchanged scalar carries no 
information about angular orientation, for each of the three two-point functions separately 
in (5.3) we can get rid of two degrees of freedom by rotational invariance (and of three by 
momentum conservation). So alltogether we get 6 — 2 — 3 = 1 degree of freedom for each 
two-point function for a total of three. This will not be the case for the graviton exchange 
discussed in the next subsection, since the tensor propagator does carry information about the 
relative orientation of the vertices it connects. Second, we would like to mention that, apart 
for the numerical coefficients, all the scalings in (5.3) could have been guessed (or computed) 
using the Feynman rules in appendix A (or the canonical formalism). 

A similar computation can be done for higher A-point functions to estimate their multi- 
vertex diagrams close to collinear limits. For example, the five-point correlation function has 
one collinear divergence (ki + k2 — > plus permutations), the six-point correlation function 
has two collinear divergences (ki + k2 + k3 — > and ki+k2 — > plus permutations) and so on 
for larger A's. It is straightforward to compute the collinear limit of an Appoint correlation 
function split into a subset with N\ external legs and its complement. We assume that N% > 3 
and N2 = N — Ni > 3 such that we can use (4.19). Introducing the notation Kn = Yli=i 

12 Notice that this expression is only valid if n B — 1 <C 3i>\/27ra, where n a is the spectral index at zeroth order 
in b. This is true in the region of parameter space we are interested in (see e.g. figure 3). 
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to leading order in b, i.e. to order b 2 , the result is 

,2x 2 Al An* At 



lim ( Kkl ■ ■ ■ KkN )<- = ^£\i% A K N » r * t E - (| " ^ ) » (I " '4) ' 

jy-h jy-h 

h<i 2 -<i h h<h-<j h ah+h hlki2 ■ ■ ■ ki h k * k * ■ ■ ■ k ih 

We would like to stress again that this result follows from systematically neglecting gravita- 
tional interactions. As we have seen in section 4.2, starting at I = 2 some terms are missed 
by this approach as compared to the exact computation including gravity. Hence only the 
leading order in a terms, i.e. l\ = 0, 1 and ^ = 0,1 are strictly correct. On the other hand, 
the subleading terms in a in (5.4) should be sufficient to estimate the combinatorics factors 
arising for large N as we will do in section 6.2. 

5.2 Collinear Limit with Graviton Exchange 

Consistency relations together with the answer (4.19) are enough to predict the behavior of 
an iV-point function in the collinear limit when the quasi on-shell internal propagator of figure 
5 corresponds to a graviton. In Appendix B, we derive the consistency relation for a collinear 
divergence mediated by a tensor mode, the main result is 

IN \' GE N!-l N 2 -l „ INt \' „ 




where GE means graviton exchange and 

Eiji m = kikjkikm sin 0i sin 6j sin 9 t sin 9 m cos(^ + 4>j - <t>i ~ <j>m) , (5.6) 

with cos#; = kj • <i/(qki), and <j>i is the azimuth angle on the plane perpendicular to q. 
Note that the derivatives are not proportional to the spectral index and they are generally 
of order one. Since we have a tensor two-point function instead of a scalar one, the GE 
contribution is suppressed 13 with respect to the SE contribution by the tensor-to-scalar ratio 
r. Also note that in (5.5) one must replace k^ = \ ^2^~ 1 k;| which introduces a dependence 
on the orientation of all momenta. Let us see how the counting of degrees of freedom in (5.5) 
differs from the case of scalar exchange. The left-hand side should depend on 3A^ — 6 real 
parameters (as usual due to translational and rotational invariance). On the right-hand side 

13 Despite this suppression, it has been argued that the collinear limit of tensor may well be the dominant 
piece in slow-roll model of inflation [8]. 
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we have one parameter for the tensor power spectrum q plus 3(Ai — 1 + A2 — 1) for the angles 
QiU<l>i,l and the norms kn. Then we have still to subtract one parameter since adding an 
arbitrary constant to all (fiij leaves the result invariant due to a cancellation in the argument 
of the cosine in (5.6). Altogether we have 3A — 6 parameters on the right-hand side as well. 
This reflects the fact that, in contrast to the scalar exchange, the graviton exchange carries 
information about the relative orientation of the split correlation functions and therefore, 
even in the collinear limit, the full dependence on the angles is preserved. 

Using (5.5) for the trispectrum leads to the particularly simple result 1 (see Appendix B 
for more details) 

% (I!**) = sin2 * sin2 *3 «* [2(fc - fc)] (^) «>' (Tliy (KlJ 

~ sin 2 d 1 sin 2 9 3 cos [2(0! - &)] ^ (K& <^ 2 3 )' {KlJ , (5.7) 
where we have use that n s — 1. 

6. Appoint Correlation Functions 

In this paper we have outlined a strategy to calculate the leading contribution to any A-point 
function in the two parameters b 1 and a S> 1 for resonant inflationary models. The main 
result is given by (4.19). We have computed more terms (4.20) but we expect additional 
gravitational contributions starting at the same order as I = 2. 

All corrections to this formula (higher orders of single- vertex and multi- vertex diagrams) 
are suppressed by extra powers of b or 1/a but might be enhanced by two effects. First, 
in some special part of k-space, a large ratio of momenta might arise, such as in collinear 
limits. This is not a serious obstacle because exactly in these regimes we can make use of 
consistency relations. Second, the number of subleading diagrams grows with A making their 
sum a potentially large contribution even if each one of them individually is very small. This 
is related to a well known problem in particle physics, where calculating scattering amplitude 
for large number of external legs is made computationally challenging by the fast growing 
number of Feynman diagrams. New techniques are being devised to handle such calculations, 
e.g. in gauge theories, and it might be interesting to investigate whether any one of them 
could be employed in the context of cosmological perturbations. 

14 This result was first obtained in [8] . There one specific direction for the collinear limit q — > was chosen 
by fixing 6\ — 63 = n/2. Here we leave 81 and O3 unspecified. Notice that (5.7) depends on six degrees of 
freedom, as it is the case for a general trispectrum. 
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In the following, using the results of the previous two sections, we estimate both of these 
effects. We find that combinatoric factors grow fast and can overcome the b and a expansion 
for sufficiently large N. We estimate the largest N, which we call iV max , for which the b and 
a expansion are trustable (summarized in table 1). For N < iV max , all correlation functions 
can be accurately computed using the results given in this paper. In the bulk of momentum 
space, i.e. for generic momenta, the result is given by (4.19). In certain localized regions of 
momentum space, which we indicate as collinear regimes, consistency relations can be used 
to compute the leading behavior. As an example we do this in detail for the trispectrum in 
subsection 6.4. 

Before proceeding, let us stress that the fact that we are not able to compute Appoint 
correlation functions above a certain iV max does not imply a failure of the general strategy 
outlined in section 2. The problem lies somewhere else. Even in the relatively simple scalar 
theory obtained after neglecting gravitational interactions, very large Appoint correlation 
functions are computationally intractable. 

6.1 Estimate of the combinatorics of single-vertex diagrams 

As we have seen in subsection 4.2, only the leading and first subleading order in a of our 
computation of single- vertex diagrams (4.19) correctly reproduces the exact computation in- 
cluding gravity. Hence our result is trustable as long as the higher orders are indeed negligible. 
Although these terms are suppressed by powers of a they are enhanced by combinatoric fac- 
tors which grow with N. In the following we estimate the smallest N, which we will call A^ ax , 
for which the I = 2 subleading terms start to dominate over the I = 0, 1 terms of (4.20). 

Using (4.20), we take the ratio of the subleading terms in a, i.e. those with I > 2, to the 
leading term plus first subleading term in a, i.e. the one with I = and 1 = 1. To simplify 
the estimate, we will assume that all fe's are of the same order, so that K = k{ ~ Nk. 
Once this is done, one finds 



A similar simplification can be done for the sum of I > 2 terms to write them all as a single 
overall amplitude times a sine up to a phase. Since this leads to cumbersome formulae, we 
prefer to simply assume that all Z > 2 terms interfere constructively, i.e. we approximate 
sin ~ cos ~ 1 in summing them up. This gives a conservative estimate of A r max , which we 
have checked is practically undistinguishable from the one obtained by correctly accounting 
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for all trigonometric functions. So our estimate goes as follows 

|1 > 2 subleading single-vertex] Sz=2 fr YliiK^-Ki^h^h ■ ■ ■ kit) 1 sin (jf~ ~ J 
I Z = 0, 1 leading single-vertex! \ _ I Y\ K cos (£k) 

~ ^= 2 ^j ill = v °i °i (6 1) 

where the binomial coefficient is given by (f) = N\/[l\(N - 1)1]. The first TV and Z for which 
this quantity becomes of order one depends on a. Both theoretical and phenomenological 
arguments guide our attention to values of a in the range 100 — 1000. Let us go through these 
arguments. 

For an inflationary background of the large-field type as in [21, 22, 23], we expect \J2e* ~ 
O(.l), while for / (remember in our conventions we have set the reduced Planck constant 
to one) we expect some subPlanckian value, i.e. / ^ 1 [46]. Since b is related to some 
non-perturbatively generated scale A, we assume b < .1. On the phenomenological side, 
we find again that b <C 1 is required by comparison of the resonant spectrum with the data 
[44, 22]. Then, the resonant bispectrum can be estimated by [18, 25, 26] / rcs = 3&v / 27ra 3 / 2 /8 ~ 
bc?l 2 . Therefore a detectable bispectra will require a of order 10 2 or larger. Putting these 
observations together we conclude that an interesting range of values is a ~ 10 2 — 10 3 . 
Demanding that (6.1) is less than 20% gives N^ a ^ for a given a. Some explicit values are 

{a = 100, iV^ = 9 }, {a = 500, N^, x = 18}, {a = 1000, = 26} . (6.2) 

6.2 Estimate of the combinatorics of multi-vertex diagrams 

The leading order contribution (4.19) only includes the single vertex diagrams depicted in 
figure 1. Each vertex comes with a factor of b and therefore the amplitude of multi- vertex 
diagrams is subleading in that parameter. Also, although it is not as obvious, the amplitude 
is subleading in a. In the following we estimate the combinatoric factors associated to these 
diagrams which could invalidate the b and a expansions. We do not directly compute any 
multi-vertex diagram, but we use the results of section 5 for the collinear limit to estimate 
them. Our final result will be an estimate of N™ x (a,b,r), i.e. the largest N for which the b 
and a expansions are not invalidated by the large number of multi- vertex diagrams. 

We start by considering only the case of scalar exchange (SE) and we look at the multi- 
vertex diagrams arising from splitting the Appoint correlation function into N\ external mo- 
menta on the one side and N2 = N — N\ on the other. For simplicity we take all momenta 
to be the same ki ~ k and adopt the notation = = Nk. We also assume that 
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N\ > 3 and N2 > 3 such that we can use (4.19). One can check that the case N\ 
N2 = 2 give a smaller contribution. We then find 
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(6.3) 



The numerical coefficient grows quickly with iV and is maximized for N2 = N± = N/2 where 
it is N 3 2 N ~ 6 (1 + N A /a 2 ) 1/2 . This is expected since our leading contribution is a single 
term in k-space while as many as iV 3 terms contributes to the collinear limits. 

An analogous estimate can also be obtained for the collinear multi- vertex diagrams with 
graviton exchange (GE). Starting with (5.5), let us focus for concreteness on the N± part, 
since the N2 will be identical. Again we assume N\,N2 > 3, since we checked that Ni = 2 
and N2 = 2 give a smaller contribution. We notice that the derivatives are with respect to 
k; • k, with i,j = 1, . . . , (iVi — 1), meaning that should be re-expressed in terms of the 
other (jVi — 1) momenta as k/Vi = — Si^ 1 k-i- The dominant term in (5.5) arises when the 
derivative acts on sin(a) inside (lZ Nl ^ giving 
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for any i, j = 1, . . . , (N12 — 1)- Contracting this with E^^i introduces a complicated angular 
dependence. We decide to give a very conservative estimate, by neglecting the cancellations 
occurring due to the various sinusoidal dependences and setting sin ~ cos ~ 0(1). This gives 
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where in the last step we have assumed for simplicity that all momenta are of the same order, 
hi ~ k. Putting things together our estimate for (5.5) is 
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where we have introduced the tensor-to-scalar ratio defined in (B.17) and the extra factor of 
two comes from the sum over polarizations. We therefore obtain the result 



The numerical coefficient is maximized for N\ = N2 = N/2 where the combinatorics is 
(N/2 - l) 2 N 3 2 N ~ e (1 + N 4 /a 2 )~ 1/2 . Comparing with (6.3), we see that the combinatorics is 
worse than for the scalar exchange by a factor (N/2 — l) 2 but there is an extra suppression 
by a factor of r which is at most a fifth but could be much smaller. We remind the reader 
that we have assumed total constructive interference between various oscillating terms and 
therefore this is a very conservative overestimate. 

6.3 iV-point Correlation Functions 

Now that we have either computed or estimated all types of contributions to a generic tree level 
iV-point function, we are ready to collect our findings into a final result. We will first discuss 
a generic iV-point function away from collinear limits, i.e. away from the localized regions 
in momentum space where a combination of momenta vanishes 15 . Thanks to consistency 
relations we will be able to return to those collinear limits later on. 

Away from collinear divergences, three quantities control the size of a given diagram: 
b, a and the combinatorics 16 . In section 4 we have obtained the leading order result in b 
and a, given in (4.19) (with coefficient given by (4.17)). Although this result was obtained 
neglecting gravity, it is a very accurate approximation to the exact result because gravitational 
interactions would give small corrections at the same order in b and a. Intuitively this 

15 Our formula also does not capture multiple squeezed or multiple collinear limits where more than one ex- 
ternal momentum (or more than one internal momentum) go to zero. We have not attempted to study this case 
thoroughly and we leave this step for future investigation. It should be noted that multi-squeezed/collinear- 
divergent terms become large only in very small regions of momentum space, which have effectively codimension 
two or more. 

16 In principle one should also include the tensor-to-scalar ratio r appearing in the contribution from graviton 
exchange. We provide a conservative estimate assuming r ~ .1, which is the largest order of magnitude allowed 
by the data. Our results are modified very little if one assumes a different fiducial value for r. 
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Single vertex 


Multi-vertex Scalar 


Multi- Vertex-tensor 


a\b 


any 


a\b 


0.1 


0.01 


a\b 


0.1 


0.01 


100 


9 


100 


12 


15 


100 


13 


16 


1000 


26 


1000 


18 


21 


1000 


18 


21 



Table 1: Estimate of the value iV max for different cases. For single vertex, iV max is defined as the 
smallest N for which the b- and a-subleading terms (I > 2 in (4.20)) add up to 20% or more of the b- 
and a-leading terms given in (4.19). For the collinear scalar and graviton exchange, -/V max is defined as 
the smallest N for which the conservative estimate of the combinatorics in (6.3) and (6.7) with a ratio 
of scales k/q ~ 1 add up to 20% or more of the b- and a-leading terms. For the graviton-exchange 
diagrams we have taken the tensor to scalar ratio r to be 0.1 in order to give a conservative estimate. 
For inflationary backgrounds of the small-field type r could be much smaller. 

is because the relevant interactions come from the scalar field sector, while the gravitational 
ones are much smaller. Therefore, as long as b <C 1 and a ^> 1 are good expansion parameters 
any iV-point function away from collinear limits is given by (4.19). 

Due to the fast growth of combinatoric factors for N > iV max , b and a cease to be good 
expansion parameters, i.e. higher orders are not negligible anymore. In the two previous 
subsections, we have estimated this growth for some b- and a-subleading diagrams. The 
smallest N for which the combinatoric factors overcome the b and a suppression are collected 
in table 1 for various diagrams and values of the parameters. For concreteness we have written 
down the smallest N for which the b- and a-subleading terms from each class of diagrams 
are 20% or more of the b- and a-leading terms 17 . The smallest N for which the sum of all 
subleading terms is 20% or more of the b- and a-leading terms is iV max . This is visualized in 
figure 3. For example for a = 100 and any b < .1, iV max = 9, while for a = 1000 and any 
b < .1, iV max = 17. To summarize, (4-19) is a very accurate approximation to all N -point 
correlation functions for 2 < N < iV max away from collinear limits. 

Let us now turn to the discussion of collinear limits. The momentum dependence of 
any iV-point function is such that divergences arise when some momenta or a combination of 
them go to zero, familiar examples being the A; -3 scaling of the spectrum of a free field in de 
Sitter or of local non-Gaussianity in the squeezed limit. In our computation, one might be 
worried that when such divergences appear in a subleading term, the latter can become very 
large and invalidate our result. Luckily for us, the squeezed limit in which a single external 
momentum goes to zero is correctly captured by the first subleading term in (4.19). On the 
other hand, collinear and multiple squeezed enhancement appears in subleading terms that 

17 Note that by leading terms we mean all the terms in (4.19). It includes the first 1/a corrections. 
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Figure 3: The plot shows N max for different regions of the log 10 (6) — log 10 a plane. As expected 
smaller b and larger a extend the validity of (4.19)-(4.17) to higher and higher N's. 

our method cannot compute. Yet again, we are lucky since anytime these singularities are 
present, the iV-point function can be rewritten in terms of lower- iV-point functions away from 
these special point using the consistency relations reviewed in appendix B. Therefore it is 
straightforward to include these collinear singularities into our result (4.19). As an example, 
we include collinear singularities for the trispectrum in the next section. 

6.4 The Trispectrum 

Even if the combinatorics are under control, one expects that the leading formula (4.19) will 
fail near collinear singularities. It would be interesting to have a better approximation of 
the iV-point correlation functions that is valid everywhere in k-space. This is in principle 
easy to do since we can use the consistency relations to find out the correct behavior, with 
coefficients, in each of these collinear limits. In practice, it is very cumbersome to write down 
all the possible terms at large N but in a given application, one will usually need only a 
specific limit which is easy to figure out from our formula in case by case situations. Here we 
will give the trispectrum including the corrections in the collinear limits (both from scalar 
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and tensor exchange). 

A very good approximation to the trispectrum can be obtained by taking our leading 
results (4.19) (with a = <j) K /f) 




-4, 



/ 2^6a 7 / 2 (2vr 2 A|) 3 , 



sin (a) 



(6.8) 



and adding the relevant momentum singularities using the various consistency relations. For 
the four point function, we have collinear divergences mediated by scalars and tensors as 
well as squeezed divergences. The latter is already included in (6.8). We also have multi- 
squeezed-limits and the case of collinear and squeezed limits taken at the same time. We will 
not include them in the present discussion as the necessary consistency relations are harder to 
derive. We leave the analysis of these terms for further work, noticing that they are relevant 
only in very localized regions and do not affect the bulk shape of the iV-point function. Using 
the consistency relation in appendix B, (4.19) for the N > 3 and (3.16) for N = 2 we find 
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(6.9) 



Note again that we have not included the double squeezed limit and that this formula is only 
valid near collinear singularities, i.e. when any one of ^1^,^1,4,^2,3 goes to zero. One should 
not take a squeezed limit of this expression, the answer in the squeezed limit is given by (6.8). 
The size of the various terms can be estimated by assuming that all momenta are of the same 
size except the vanishing ones. One finds 

squeezed 4 cos(a) k 



a sin(a) kj ' 
24\/27r6 sin(a) fc 3 



a leading 
collinear — scalar 

a leading a 5 / 2 

collinear — tensor 3r sin 2 9\ sin 2 0% cos(2(^>i + ^3)) k 3 

P~ 

ft 12 



ju3 
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(6.10) 
(6.11) 
(6.12) 



a leading 2 v / 27r6a 7 / 2 sin(a) 

where kj and k\2 are g° m g to zero. We see that all these terms are very small in the equilateral 
limit, since the combinatorics are not yet important here, but they are enhanced by a large 
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ratio of scales near collinear singularities. For CMB observations this ratio can be as large as 
k/kj ~ k/k 12 ~ 10 3 . We can estimate the size of the trispectrum by evaluating it at different 
points, say the equilateral, collinear and squeezed points, 



ki ^^4) (equilateral 
TO 



1 /v k4/ lequilateral /„ 1 Q \ 

*NL = f^Tje , (6.13) 
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9NL = — A ,„ |2|p |4 (6.15) 



where |i?fc| 2 is the power spectrum define in Eq. (3.16). Note that our definition of ijvz, differs 
from the one in [7] by a factor of 8 while tnl and qnl agrees with the literature [50] when 
taking the appropriate limits. Evaluating our trispectrum in the appropriate limits using 
Eq. (6.8) almost everywhere and Eq. (6.9) in the collinear limits we get 

sin a 

tNL = — g , (6.16) 

9 

t nl = — r sin 2 #1 sin 2 (9 3 cos [2(0i + 3 )1 + 18?r6 2 a sin 2 (a) , (6.17) 
16 

75 

9NL = -— \/2^6a 5 / 2 cosa . (6.18) 

lUo 



Note that we have not specified exactly how one takes the limit k±2 — > to get tnl which is 
why there is still an angular dependence. Current bounds on the trispectrum can be found 
in [51]. Since our signal is oscillatory, probably a dedicated analysis is requires to extract it 
from the data [45, 44]. Notice that in writing the above formulae we are omitting an arbitrary 
phase in sin(a) and cos(a), which is not fixed by the model. 



7. Discussion 



In this paper we have computed for the first time the iV-point primordial curvature correlation 
functions from inflation up to N of order ten or more at tree level. We did this in resonant 
inflationary models [18] with the potential (3.1), which find a natural UV completion in 
models based on axion monodromy [20] constructed in string theory [21, 22, 23]. Our main 
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results is 



N<N, 



max 



( _ ) iV£^ a 2iV-9/2 A 2iV-2 (27r 2 ) 



(2tt) 3 5 3 5>i A N B N , 




(7.1) 




This formula captures the leading behavior of the N-point correlation away from collinear 
singularities for N < iV max , where N max is of order ten or larger depending on the value of 
the parameters. Corrections to this formula are less than 20% for any iV < iV max (see figure 
3 for the list of iV max ). Let us review the most important ingredients of our calculation. 

First, we have noticed that the main interactions in our inflationary model come from the 
scalar field sector and can in principle be made parametrically large (by taking a large a). This 
justifies neglecting gravity, which drastically simplifies the computation. We have checked that 
our result is an excellent approximation to the exact result including gravity, reproducing it 
for N = 3 at leading and first subleading order in a. The success of this decoupling strategy 
was already noticed several times in the literature and ours is an additional motivation to 
formalize it for a generic case. 

Even after the problem has been simplified to computing correlation functions for a 
relatively simple scalar field theory, some further obstacles are encountered. The first is the 
presence of collinear limits (or multiple squeezed/collinear limits). The issue is that even very 
small contributions to a correlation function, which are potentially sensitive to gravitational 
interactions, could in principle be boosted by momentum enhancement arising in some par- 
ticular limits. Since our strategy is to decouple and consequently neglect gravity, these terms 
seem hard capture. Fortunately, consistency relations exist that relate squeezed/collinear 
limits to correlation functions with less external legs, the latter evaluated away from these 
points. The second obstacle is the fast growth with N of combinatoric coefficients. This is 
a well known problem in field theory and in fact a very active field of research. We have 
not been able to overcome this obstacle which therefore sets the limits of validity of our 
computation defining an -/V max above which our result (7.1) is not trustable. 

Since our approach has been partially inspired by the effective field theory of perturba- 
tions around an inflationary background proposed in [17] (and recently generalized to mul- 
tiple fields in [52]) let us comment on the connection. The authors of [17] start from the 
most general Lagrangian for perturbations around a homogenous time dependent inflation- 
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ary background. In order to recover time diffeomorphism invariance one can perform the 
Stiickelberg trick introducing a new field tt, which non-linear ly realizes the broken symmetry. 
This formulation stresses that just the effective interactions amongst the perturbations are 
relevant for actual observables so that one is able to differ a precise construction of what 
generates the inflationary background. When the underlying theory is specified as in the res- 
onant model considered in this paper, then it is straightforward to switch from the standard 
description of perturbations to the tt gauge. Independently of one's favorite variables, the 
important point is that there is a decoupling limit in which the metric can be approximated 
as unperturbed. For observables not dominated by gravitational interactions this limit can 
be used (as we did) in order to drastically simplify computations. To clarify the equivalence 
between the description of [17] and the one used in the present paper, let us briefly sketch 
how our computation translates into the tt variable. Since our model has a canonically kinetic 
term and no higher derivatives 18 we can set M2 = M3 = in [17]. One can then follow the 
derivation in [53] and solve the constraints of the most general gauge-unfixed Lagrangian for 
g^ u and tt using the ADM formalism. At linear order one finds 19 

5N = eHir, (7.2) 
diN' 1 = -eH-k + e 2 H 2 7r . (7.3) 

No slow-roll approximation was used to obtain this result. In the action there are many 
couplings between N, N l , the metric perturbations 7^ and tt. Substituting the value of the 
lapse and shift found above into the action, one can systematically check that the kinetic term 
(dir) 2 dominates over everything else at the quadratic level. This is because in resonant models 
e <C 1, so that, at least at the level of linear constraints, decoupling works as in standard slow- 
roll inflation i.e. one can simply study the decoupled tt Lagrangian approximating the metric 
as unperturbed. Then the crucial point is to keep higher derivatives of the time dependent 
coefficients (in our case only H and its derivatives) when expanding in small tt. For example 
for the bispectrum the relevant terms is the one with the largest number of derivatives, i.e. 
— H^tt 3 /6. To leading order in a this is the same as V"'5(ft 3 /6 given the relation Sep = —Tr<j) 
valid at linear order. 

18 The simultaneous presence of modulations in the potential and higher derivatives has been studied in [56] 
using the P(X, <f>) formalism. 

19 This result was derived in [2] both in comoving and spatially flat gauge. In the presence of higher powers 
of single derivatives the right hand side of (7.3) is simply multiplied by cj 2 , while the constraint for N is 
unchanged. 
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A systematic term by term inspection of the decoupling limit can also be done using the 
54> action derived in [2]. Every terms in Eq. (3.8) of [2] except the V" term are gravitational 
in origin. By inspection, they are all subleading with respect to V'" for resonant inflationary 
models. 

To conclude, we mention some interesting directions for future research: 

• Determine a general formulation of the decoupling limit and provide some formal proof 
of its validity. 

• So far all examples seem to confirm the validity of the strategy outlined in section 
2. Determine whether there are counterexamples. Can deviations from Bunch-Davies 
[54, 55, 56] be seen as single field models with parametrically large non-Gaussianity? 
Can we have a simple quantitive test of whether a theory decouples or not? Models 
with features in the potential are also single field models with parametrically large 
non-Gaussianity, can the decoupling limit be applied to them? 

• Investigate in detail consistency relations for multiple independent squeezed/collinear 
limits. 

• Use our results to perform a detailed phenomenological analysis of the trispectrum (and 
higher point functions) for resonant models. 

• Investigate the mechanism of resonant enhancement for isocurvature perturbations. 
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A. Feynman rules 

In this appendix, following [47] we derive the Feynman rules for the theory for the perturba- 
tions that we study in the main text. The expectation value of an operator Q at time t > ti 
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that starts in the \in) state subject to a time dependent Hamiltonian H{t) is 



(Q(t)) = (m| 



Texp ^i J dt'H(t')^ Q Texp^-ijf dt'H(t' 



We can also calculate correlators using the path integral 

rt 



Z[J+, J_] = y V4> + V<jTexp i jf dt' y d 3 x (£[#+] - £[0~] + J+</>+ + J. 



(A.l) 



(A.2) 



Correlation functions are then obtained by varying Z with respect to the sources. Differenti- 
ating with respect to J+ brings down a time ordered <f> while differentiating with respect to 
J_ brings an anti-time ordered <p. Because of the doubling of fields, the Green's functions can 
be arranged in matrix with four entries containing all the different time ordering G^(x, y). 

( G++(x,y) G+~(x,y)\ 
\G-+(x,y) G+ + (x,y) ' 



(A.3) 



G-+(x,y)=i(<f>(x)<t>(y)) 
G ++ (x,y) = i(T ( f>(x)^(y)) 



G + -(x,y) = i(<P(y)^(x)) , 
G (x,y) = i (T(/)(x)<p(y)) . 



(A.4) 



They obey the identity G ++ (x,y) + G (x,y) = G + (x,y) + G + (x,y). We will use a 
different basis where 



+ 



4>c 



and where the Green's function are given by 



iG c G R 
G A 



R 



G ++ G H 
G-+ G~ 



R 1 



(A.5) 



(A.6) 



for some rotation R, with 



Gc(x,y) 



{G- + (x,y) + G + -(x,y)) 



G R (x,y) = G ++ (x,y)-G + -(x,y) , 
G A (x,y) = G ++ (x,y)-G- + (x,y) . 



(A.7) 



We can build the various propagators from the modes solutions given in (3.15). It is useful 
and customary to Fourier transform the space components while leaving the explicit time 
dependence. 



Gc(fc,ri,r 2 ) = Re (H«(-kTi)H«*(-k7a) 

4a(n)a(r 2 ) V 

G a (k,T l3 ^) = - e(n - r 2 )Im (h« (-krOHW* (-kr 2 



(A.8) 
(A.9) 
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with a(r) ~ —1/(Ht) to first order in slow-roll expansion and v = 3/2 + 2e + 8. We will 
denote the field <pc by a solid line and the field (ft a by a dashed line. To find the Feynman 
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Figure 4: Propagators and some vertices (AT = 3 and N = 4) for this theory. Note that G J 4(fc, n, T2) = 
Gfi(k, T2, ti). We use the shorthand notation ^12 = <5(ti — T2). When a two-point function is attached 
to a vertex, the corresponding time must be integrated over. Internal spatial momenta must also be 
integrated over J d 3 p/(27r) 3 . As usual, one must also divide by the symmetry factor for the diagram. 
At A" = 5, there are 3 diagrams, etc. 



rules we rewrites the potential in terms of 4>c an d 4>a where for simplicity we denote Sep = (/), 

E^a^-w-n-E^Ef^^-Vi^-H*) , (A.10) 



N=3 N=3 k=0 



where is given in (4.8). In our Feynman rules, the coefficient of a given vertex is the 

coupling appearing in front of a term 

rh N ~ k (h k 

A N-Kk {N _ k)lkV (A.llj 
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which is simply (for k running from to N) 

4&_ M = (1 ~ 2 l" )fc) ■ (A.12) 

Note that this vanishes for any even k. See figure 4 for the propagators and the first few 
explicit vertices. 

B. Consistency Relations for Squeezed and Collinear Limits 

In [2] it was realized that, in single field inflationary models, interesting consistency relations 
arise in the limit where one of the momenta of the three-point correlation function is taken 
to zero (squeezed or soft limit). As it was later stressed in [27], the key property is that there 
exists a single clock on which all correlation functions depend which is valid for all single 
field models with no assumptions about slow-roll dynamics. The original derivation was done 
for the bispectrum but the generalization to any iV-point functions is straightforward [6, 29]. 
Here we re-derive it in x-space following [53]. 

With very analogous reasoning one can also derive a consistency relation for collinear 
limits, i.e. when a subset of the momenta of the external legs (and its complement) sum up to 
zero. This relation has also been used in the literature before [8, 30, 57]. We use the notation 
that prime in () means that the factor (27r) 3 <5 3 (K) is absent. 

Squeezed (or soft) limit 

Correlation functions in single field inflationary models depend on a single clock and are 
evaluated at horizon crossing. Using the inflaton as a clock, all iV-point functions should 
be evaluated at <j> = where 4>* corresponds to some pivot scale fc* relevant to a specific 
observation. In a squeezed or collinear limit of a given mode q — > 0, the associated curvature 
perturbation TZ B , where B stands for background, is well outside the horizon, hence it becomes 
classical and its constant value is determined by the square root of the spectrum. The trick to 
evaluate an iV-point function in this limit is to perturb the background by this long wavelength 
mode TZ B as in (3.8). Assuming a small perturbation 1Z B , by Taylor expanding one finds 
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A shift in 1Z B can be traded for a shift in time, or 0* or fc* using dlZ = Hdt = d^^j^flt* 
tflogfc*. If we decide to use cp* to parametrize time, we find 
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(B.2) 



The squeezed limit is then obtained by multiplying the last expression by some IZk and taking 
the expectation value 
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lim/ftqlJ^ \ = (2ir) 3 5 3 (K)V2e;\n q 
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(B.3) 



i=i I \i=i I 

where it is understood that everything is evaluated at horizon crossing. 

Collinear Limit, Scalar and Graviton Exchange 

A similar consistency relation can be derived in the collinear limit where an intermediate 
momenta goes to zero, see figure 5. The intermediate state could be a scalar or a tensor 
and we will discuss both. The first correlation function for which such a limit exists is the 
trispectrum [8, 30]. It is straightforward to find a relation for any TV-point function. Consider 






Figure 5: The behavior of an TV-point function in the limit where an intermediate scalar or graviton 
propagator has vanishing momenta can be obtained using the consistency relations. 
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the case in which a proper subset of the iV external momenta, say 1, • • • , r, sum up to q and 
we take the limit q — )■ 0. The other momenta r + 1, - ■ ■ N sum up to zero as well by momentum 
conservation, such that Y2l=i ki = Y2jLr+l = Then diagrams like the ones in figure 5 
have an internal propagator going on shell and hence can become very large. This quasi-on- 
shell propagator that mediates a long-range interaction (q — > 0) can be either a scalar mode or 
a graviton mode. In the following we will treat each case separately. Let us start considering 
the scalar-exchange (SE) case. We can separately evaluate the correlation function of the 
first and second subset of momenta in the background of some long wavelength scalar mode 
Kf. For small K B we can Taylor expand and use dlZ = d(f>/y2e as in (B.2). The consistency 
relation in the collinear limit is then obtained by multiplying these two expressions and taking 
the expectation value. The result is 




(B.4) 



Since \1Z q \ 2 oc q~ s , this contribution is divergent in the strict collinear limit at tree level. 

Let us now turn the case in which a long wavelength graviton is exchanged (GE). This 
situation was first analyzed in [8] for the four-point function. We follow their derivation 
generalizing it to any iV-point function and we will work in momentum space. Analogously 
to the case of a scalar mode, we consider a TV-point correlation function in a background 
perturbed by a long wavelength classical graviton, which appears in the metric at linear order 
as 

9ab = a 2 (5 ab + 7 ab ) , (B.5) 

where 7 is transverse and traceless. First we want to translate derivatives with respect to 7 
into derivatives with respect to fcj. Using the perturbed metric (B.5), we find 

k 2 -> k a k a + k a lab k b = k a k a - k al ab k b (B.6) 

where a,b = 1,2,3 and the sum over repeated indices is implicit. Scalar functions, like the 
iV-point scalar correlation functions, due to rotational invariance can depend on (spatial) 
momenta only via combinations like k a g ab k b . The correlation function (7^ki ' " • T^U N Y is a 
function of {ki, . . . , kjv_i} and so we can use the chain rule 



Then, Taylor expanding we find 
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For a massless spin-two field like the graviton, the two propagating degrees of freedom can 
be put together into a causal quantum field as 



7ab (k)= £ [ 7 Ut)ble: b (k) + rt(t)b%e: b (-k) 

s=x,+ 

(2^) 3 <W5 3 (k-kO, 



(B.9) 



where e* fe (k) are the two polarization tensors satisfying e^ a (k) = 0, k a e s ab ik) = as well as the 
completeness relation e^ b (k)e^ b (k) = 2S ss i. Now, as we did in the scalar case, we can multiply 
an Ni and A^-point correlation functions in the background of some long wavelength graviton 
and average over the background. The result written in terms of polarization tensors is 
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In order to perform the sum over polarization, we choose an orthonormal basis {e 1 , e 2 , q/q}. 
The polarization tensors can be written as 



eW - e 2 e? e x 



e^e? + e 2 ei . 



Using spherical coordinates, each vector can be expressed as 

kj = fcj(sin ftj cos fa, sin ftj sin fa, cos ft 



(B.ll) 



(B.12) 



where cosftj = kj • c[/(qki), and fa is the azimuth angle on the plane perpendicular to q. A 
straightforward computation then shows that 

e+ b (q)fcf fcj = kikj sin ftj sin Oj cos(fa + , 

e* 6 (q)fc"fc- = kikj sin ft> sin Oj sm(fa + fa) . (B.13) 
All together, the polarization sum is 

Eijlm = € lb( c l) e cd( c l)Kk b jkik'? n 

s=+,x 

= kikjkik m sin 9i sin Oj sin ftj sin m cos(fa + 0j — fa — 4> m ) . (B.14) 
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Hence we find the general expression 
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(B.15) 

Note that it is understood that the correlation function () have their last momenta eliminated 
using translation invariance (e.g fcjvi = I Sf^ 1 ^!)- This means that even if (7t N \ is only 
a function of the norm of each vectors (as we have in (4.19)), (lZ N ) depends on angles and 
the derivatives get quickly very complicated. The special case is N = 4 with N\ = N 2 = 2 
where the formula simplifies to 
1 a \ iGE 



a 



lim ( n^ k * ) = hi\ 2k i k2 s sin2 ^ sin2 #3 cos [2(0i - 3 )] ^ (K^K-^)' ^ <^k 3 ^-k 3 )' . 

(B.16) 



U=l 



Using the standard definitions of scalar spectral index and tensor-to-scalar ratio 
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for some pivot scale k*, we find 
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lb 

where we have use n s « 1 in the last line. 

x-space Derivation of the Consistency Relations 

In the following we provide a derivation of the consistency relation (B.3) and (B.4) starting 
from rr-space. We follow the derivation presented in 20 [53], generalizing it to any iV-point 
functions. Let us start considering the connected TV-point function (lZ N ^ in coordinate space 
in a quasi de Sitter background slightly perturbed by a small long-wavelength classical mode 
TZb- By Taylor expanding around TZb = 0, we find 

d 



n 



N\ 



^ 7V ) (xi,...,x 7V ) =ft B (x) 



(xi,...,Xtv) 



(B.19) 



n B =o 



Obtained in collaboration with P. Creminelli. 
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Since the background mode TZb is very slowly changing (in the sense that its momentum q 
is much smaller than the momenta associated with TZ(xi) for i = 1, . . . , N), we can evaluate 
Kb at any one of the points (xi, . . . , xjv) or a linear combination thereof with coefficients 
of order one. For concreteness and later convenience we choose to evaluate it at xjy- The 
derivative with respect to TZb can be changed into logarithmic derivatives with respect to the 
coordinates by noticing that an additive shift in 7Z is equivalent to a multiplicative shift in 
all coordinates (see e.g. (3.8)) 

d 



d!Z B 



N 

E x ^ v 



j ■ 



(B.20) 



where Vj is the gradient with respect to Xj . Performing a Fourier transform on both sides of 
(B.19) from (xi, . . . , x^v) to (ki, . . . , kjv), we get 
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E x ^ V i<^>7e fl (xi,---,**) 

i=i 



K B =o 



Because of the translational invariance of the background, the correlation functions do not 
depend on all of the 3N real parameters (xi, . . . , xjv), but just on 3(iV— 1) of them. One simple 
way to impose this is to change variables to {xj, x^r} with Xj = Xj — x^r for i = 1, . . . , iV — 1. 
Then the correlation functions may be chosen to depend just on {xj} and not on x^v- The 
fact that this coordinate is the same where we chose to evaluate TZb simplifies considerably 
the derivation. Changing coordinates and performing the simple d^xjy integral, we find 
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where k^ = 1 kj. These identical 3(iV — 1) integrals can be performed using 





dxe- txk x^f(x) 
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The result is 



in 



>n T - 



(K N ) (k 1 ,...,k N ) = -lZB(k t 
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(B.22) 
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where now Vj stands for a divergence with respect to k,-. If we now multiply both sides by 
IZq and average over it, we find 



^(^ 7V »(q,k 1 



lim (KaR. N ) 



(B.23) 
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3=1 



quod erat demonstrandum. The collinear relation is also simple to obtain with the correct 
overall 5 function. 

We can further massage the results by noting that X^^o. 1 Hj'^j ^ s ^ e dilatation operator. 
Using the fact that the TV-point function is of degree —3(N — 1) and defining a single pivot 
scale k* with which we can form dimensionless ratios, we can further simplify the results. 



We can write (with angular dependence suppressed) (lZ N y = f(ki)g J where / is a 
function of degree — 3(N — 1) and g is a function of degree zero. This means that h is a 
linear function such that X^jj 1 • X7jh(ki) = h(ki). With this functional form for (TZ N ^, 
the previous expression simplifies to 

/ N \ 



(B.24) 



Ml(KJI n ) = (2tt) 3 5 3 U + X> 1 P q 
which agrees exactly with (B.3). 
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